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Abstract 

A sign of topological order in a gapped one-dimensional quantum chain is the existence of edge 
zero modes. These occur in the Z2-invariant Ising/Majorana chain, where they can be understood 
using free-fermion techniques. Here I discuss their presence in spin chains with Z„ symmetry, 
and prove that for appropriate couplings they are exact, even in this strongly interacting system. 
These modes are naturally expressed in terms of parafermions, generalizations of fermions to the 
Z„ case. I show that parafermionic edge zero modes do not occur in the usual ferromagnetic and 
antiferromagnetic cases, but rather only when the interactions are chiral, so that spatial-parity 
and time-reversal symmetries are broken. 

1 Introduction 

One of the fundamental results of statistical mechanics is Onsager's computation of the exact free 
energy of the two-dimensional Ising model [lj. This provided a canonical example of a non-trivial 
phase transition between an ordered phase, where the Z2 spin-flip symmetry is spontaneously broken, 
and a disordered phase. This computation and many others in the Ising model are greatly simplified 
by Kaufman's rewriting of transfer matrix [2] in terms of free fermions. These results apply to different 
couplings in the two directions, including the limit where the Id quantum Ising Hamiltonian can be 
extracted from the transfer matrix. In this limit, Kaufman's mapping amounts to the Jordan- Wigner 
transformation of spins to fermions [3J. However, because this mapping is non-local, computing many 
other quantities of interest requires considerably more effort 0]. A gargantuan amount of study 
therefore has been devoted to this model in the many years since these computations, and as a 
consequence, many aspects are now very well understood. 

Nonetheless, some fundamental aspects of the Ising model have only been understood well in 
recent years. In particular, Kitaev gave a dramatically different way of understanding the phases and 
the physical implications of the model |5.j in its fermionic description. The point is that a "physical" 
quantity is typically measurable locally, so the fermions in the Id quantum Ising chain are unphysical 
from the point of view of the spin system; the mapping to fermions is a very useful but purely 
mathematical trick. However, since electrons can be trapped in one dimension (for example in a 
"quantum wire"), a relevant question is: What is the physics of the quantum Ising chain when the 
fermions are the underlying physical degrees of freedom? 

Reinterpreting the quantum Ising Hamiltonian in fermionic language is very interesting but not 
very difficult. Since there is no U(l) symmetry in the Ising chain, the number of fermions is only 
conserved mod 2; in fermionic language this means there is a Cooper pairing interaction. Working 
out the fermion dispersion relation shows the unusual feature that there is only one fermi point: as 
opposed to typical models of fermions on the lattice, there is no doubling. In currently fashionable 



parlance, the fermion is Maj'oranajj 

Less simple, however, is giving a fermionic interpretation of the phases. The order parameter is 
simply the magnetization, local in the spin variables, but non-local in the fermions. In the ordered 
phase, there are two degenerate ground states, so that the Z2 spin-flip symmetry exchanging them 
is spontaneously broken. Kitaev's remarkable observation is that this Ising spin order becomes in 
fermionic language an example of topological order. A signal of topological order is the appearance 
of degenerate ground states without a local order parameter being spontaneously broken. The Ising 
chain rewritten in terms of Majorana fermions fits the bill; the local order parameter in terms of spins 
becomes non-local. 

Taking open boundary conditions is an excellent way of probing the situation. Gapless or zero- 
energy edge modes often (although not always) occur in a phase with topological order. In the 
Ising/Majorana chain, the Z2 conserved quantity is (— 1) , the fermion number mod 2. One of the two 
ground states thus can be taken to have (—1) F = 1, the other (—1)^ = —1. Any operator \E' mapping 
one such ground state to the other must necessarily be fermionic, i.e. obey (—l) F ^f = — 1)^. A 
fermionic zero mode, however, does more than just map between the ground states: it commutes with 
the Hamiltonian, so that the entire spectrum in the two sectors (— 1) F = ±1 must be the same. Since 
the phase with spin/topological order is gapped, such a zero mode can be localized at the edge of the 
system. For the Ising/Majorana chain, there are two such exact fermionic edge zero modes [5]. These 
undoubled (or "unpaired") Majorana edge zero modes are a defining characteristic of topological 
order in this system. 

This observation has attracted a great deal of attention because in the gapped topological phase, 
the two ground states can be used to make a "qubit", a quantum two-state system robust against 
decoherence. It is robust because the gap makes transitions to any excited states essentially impos- 
sible, while transitions between the two ground states can only be caused by processes that add or 
remove an electron from the system. Acting with both zero modes does not change (—1) F , but since 
the modes are localized at opposite ends of the system, no local noise can cause this change. The 
non-locality of the order parameter in the fermionic basis is thus a feature, resulting in this protection 
against decoherence. Various proposals have been made to realize the Ising/Majorana physics in real 
fermion systems, and there are reports of observation [6]. 

A natural next step is to understand the effects of including interactions between the fermions. In 
this one-dimensional system or a collection of coupled such systems, the answer is known. As long as 
the number of chains is non-zero mod 8, the only way to go from topological order to another phase 
is to close the gap, i.e. tune the parameters through a quantum critical point [7j. The topological 
order that remains is essentially the same as that of the systems without interactions. Finding 
topological order in Id systems that are "far" from free fermions (i.e. those that can not be obtained 
by deformation) is a much more complicated question. A general classification scheme exists [U |9l [10] , 
but connecting this to simple models is non-trivial. 

It thus seems a good idea to go back to spin systems to understand if any simple-to-describe 
models also can be reinterpreted in terms of topological order. A natural set of candidates are clock 
models, where the two-state Ising spin variable and a Z2 symmetry of Ising are generalized to having 
n states and a Z n symmetry. Not only do these models exhibit continuous phase transitions between 
order and disorder [TT], but they exhibit an even richer spectrum of behavior in the presences of 
chiral interactions, where spatial-parity and time- reversal symmetries are broken \12 \ W& \ PH]. Even 
more enticingly for present purposes, the clock models can be rewritten in terms of parafermions, Z n 
generalizations of fermions [15] . 

1 Surprisingly, the word "Ising" does not appear in Kitaev's paper (except inside "Surprisingly") nor in the excellent 
review article [6] ; often this chain is now referred to as the "Kitaev chain" . In an effort to restore interspecies harmony, 
I will usually call it the "Ising/Majorana" chain when written in terms of fermionic variables. 
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The purpose of this paper is to describe how edge zero modes can occur in spin chains with Z n 
symmetry. These modes are described naturally in terms of the parafermions. For certain choices of 
couplings in these systems, I prove that exact edge zero modes exist, analogous to the Ising/Majorana 
chain. This proof applies even to spatially varying couplings, i.e. the "random" case. The correspond- 
ing phases thus are presumably topologically ordered. 

One result of the analysis here is that the most widely-studied examples of clock models, those 
with ferromagnetic or antiferromagnetic interactions, do not have edge zero modes (at least by this 
construction). Rather, edge zero modes generalizing the Majorana ones only occur if spatial-parity 
and time-reversal symmetries are broken, i.e. the model is chiral. This may seem peculiar, but in the 
classification schemes for topological order in free-fermion systems in arbitrary dimensions, discrete 
symmetries under time reversal symmetry charge conjugation play central roles |16|. I17j. Breaking or 
including these discrete symmetries can indeed change the type of topological order, or eliminate it 
altogether. The same principle seems to be applicable to this strongly interacting system. 

As will be described below, chiral interactions in the clock models result simply from including 
phases in the nearest-neighbor interaction terms; in fact, one can continuously interpolate between 
ferromagnetic and antiferromagnetic interactions by changing these phases. The most robust type 
of interaction for topological order seems to be essentially halfway in-between ferromagnetic and 
antiferromagnetic. In this case the model is symmetric, having a discrete "charge-conjugation" type 
operation that takes the Hamiltonian H —> —H while leaving the spectrum invariant, just like the 
Ising/Majorana chain does. 

This analysis here is considerably more involved than in the Ising/Majorana case. The reason 
is that the parafermions are not free in any sense; one can for example not use Wick's theorem to 
simplify the computation of correlators. Thus while they are exceptionally useful for deriving the zero 
modes here, they do not automatically make computations as easy in general. In fact, for the nearest- 
neighbor interactions considered, the chiral clock chain is only integrable in a two-parameter subspace 
of the full model [T5] . This subspace does not even include the ferromagnetic or antiferromagnetic 
cases, except at the critical points [20]. Rewriting the model in terms of parafermions does however 
shed light on the integrable cases, in fact giving a simple way of characterizing the integrable cases; 

1 will return to this in a separate publication [21] . 

The outline of the paper is as follows. In section [2j I review some basic facts about the Ising 
model, including the presence of fermionic edge zero modes in the ordered phase. The Z n -invariant 
spin chains and their parafermionic descriptions are introduced in section |3j The topic of section [4] 
is the parafermionic edge zero modes, and developing an iterative procedure to find them. They are 
easily found in an extreme case of the couplings, but even at next order they only occur for chiral 
interactions. In section [5j I prove that for appropriate parameters, the edge zero mode remains exact 
to all orders. Finally, in section [6] I give the conclusions and some ideas for further work. A proof 
that the number of zero modes (not necessarily edge modes) increases exponentially with the size of 
the system is given in an appendix. 

2 Zero and shift modes in the Ising/Majorana chain 

To begin, I review the quantum Ising chain, following [3j, and how to detect the topological order 
there, following [5]. 

2.1 The Hamiltonian and the fermions 

The Hilbert space for the quantum Ising chain consists of a two-state quantum system, i.e. a "spin" 
with spin 1/2, at each of the L sites. The Hamiltonian is comprised of two types of terms, those that 
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flip a spin at a given site, and those that give an interaction energy to adjacent spins. Precisely, with 
open boundary conditions the Hamiltonian is 

#IM = -/X>J-jf>K +1 (1) 

3=1 3=1 

where the Pauli matrices o~f act non-trivially at each site i, and / and J are non-negative and real. 
This Hamiltonian has a Z2 symmetry under flipping all the spins. For reasons to be apparent shortly, 
it is natural to name this spin- flip operator (— 1) , so that 

L 

(-lf = l[a]. (2) 
3=1 

This operator indeed squares to 1, and it is easy to check that [(— 1)^, Him] = 0. Another interesting 
property of the Ising Hamiltonian is that its spectrum is invariant under sending Him —Him- This 
follows from the canonical transformation o~j — > —a? on all sites, and cr| fc — > — <r| fc on every other 
site, which preserves the algebra of o~ x and a z . 

The physics of the Ising chain is well understood. There are two phases, with a critical point at 
f = J. The phase for / < J is ordered, i.e. the two-point function {aja k ) in the ground state goes to 
a constant value in the limit \j — k\ large. This is obvious in the extreme case / = 0, where there are 
two ground states: all spins up, and all spins down. In the disordered phase f > J, there is a unique 
ground state. This is easily seen in the extreme case J = 0, where the sites are independent so that 
the ground state is the eigenstate of crj with eigenvalue 1 on each site. 

A classic result of statistical mechanics is that the quantum Ising chain can be mapped onto a 
model of free fermions [3]. The fermionic operators are given by a non-local combination of the spin 
operators known as a Jordan- Wigner transformation. This is possible because of the duality of the 
model; the fermions are defined by multiplying the spin operator by its dual, the disorder operator. 
Duality is a non-local transformation, so the disorder operator is necessarily non-local in terms of the 
spin operators. Thus the fermion is non-local in terms of the spins, and vice-versa. Precisely, at each 
site of the chain, there are two "Majorana" fermion operators defined by 




\fe=l / 



Because of the non-local "strings" attached, operators at different points no longer commute. Instead, 
they anticommute: 

{aj,a k } = {bj, b k } = 2S jk , {aj, b k } = (5) 

for all j and k. Note also that that aj and bj are hermitian and that each squares to 1. It can be 
useful to define the complex fermions cj- = aj + ibj obeying the usual anticommutation relations. 
The open-chain Hamiltonian ([!]) rewritten in terms of these fermionic variables is 

L L-l 

Hm = if^ajbj + iJ ftjQj+i • (6) 
3=1 3=1 

The Ising Hamiltonian is easy to rewrite in terms of the complex fermions. One finds it includes 
"Cooper-pairing" terms involving CjCj+i and cjc- +1 . These terms do not conserve fermion number 



4 



generated by F = ■ CjCj, but only fermion number mod 2. The corresponding symmetry generator 
(—1) F is the product of all the Majorana fermion operators: 



L 




It commutes with any product of an even number of fermion operators, while it anticommutes with 
any product of an odd number. For this reason, it measures the number of fermion operators mod 
2, hence the name. This symmetry operator is also sometimes called "fermionic parity", or simply 
"parity" for short; the latter I avoid to prevent confusion with spatial parity symmetry. Since Him 
is comprised of fermion bilinears it indeed commutes with (—1) F , and so both operators can be 
simultaneously diagonalized. The states therefore can be divided into two sectors, even and odd 
under this Z2 symmetry. 

2.2 Edge zero modes 

The fermionic formulation makes it simple to find the zero-energy edge modes characteristic of topo- 
logical order. A fermionic zero mode \& is an operator that 

• commutes with the Hamiltonian: [H, ^f] = 

• anticommutes with (—1)^: {( — 1)^, H} = 

• has finite "normalization" even in the L — > 00 limit: = 1. 

The second property guarantees that \& maps the even sector to the odd, and the other two then 
require that the spectrum in the even and odd sectors be identical. An edge zero mode has the 
additional property that it is localized near the edge. This means that the dependence of matrix 
elements of $ on the fermions in the states a distance I from the edge must be exponentially small 
in I. 

It is easy to see that an edge zero mode exists in the extreme case / = of the Ising chain, where 
the flip terms vanish. In this case, the fermions a± and b^ do not appear in the Hamiltonian with 
open boundary conditions. They therefore commute with all the bilinears in the Hamiltonian, and so 



Since a\ and hi, anticommute with (— 1) and square to 1, each satisfies the other two properties as 
well. Thus each is a zero mode, and obviously an edge zero mode as well. 

For the Ising model, there are various formal arguments showing that the edge zero modes persist 
throughout the ordered phase, and so the spin-ordered phase is the same as the topologically ordered 
one [5]. However, as also noted there, because the fermions are free, one simply can find the edge zero 
modes exactly. The easiest way to find them is to utilize an iteration procedure. Commuting H\m 
with a\ gives —2ifb±. This result can be written as a commutator by using 2bi = [6102,02]; Thus a 
combination of a\ and 02 commutes with the first two terms in the Hamiltonian: 



[H IM (f = 0),oi] = [H m (f = 0), b L } = 0. 



ifaibi+i,Jbia2, «i + —02 = 0. 



This process can be iterated by using 



[Him, aj] 
[Him, bj] 



-2ifbj + 2iJbj-! 
2ifa,j — 2iJcij + i , 



(8) 
(9) 
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where by convention bo = d£+i = 0. Thus the (almost) zero modes localized around the left and 
right edges for / < J are given by 



*left = ai + ja 2 + (~\ 03 + . . . , ^ r i g ht = b L + ^6 L _i + b L „ 2 + ■■■ (10) 

These do not quite commute with Him, but instead: 

[H JM ,% e {t]=f(j\ b L , [Fim, bright] =/(j) °1. 

The coefficient is exponentially small in L, so in the L — > 00 limit, each commutes with Him- Each 
remains normalizable for / < J, so throughout the ordered phase, each is an edge zero mode. Thus 
Ising spin order indeed translates in the fermionic language to topological order. 

An important fact is that the edge zero mode survives in the presence of couplings varying over 
space. Namely, for 

L L-l 

H = i Yl fj a 3 b 3 + * Yl J > h '"> • 1 • ( n ) 

3=1 3=1 

the left edge mode is modified to 

*left =0.1+ —a 2 + -r-ra,3 + ... , 
J\ J\J 2 

and analogously for the right. As L — > 00, these commute with the more general Hamiltonian. They 
remain normalizable throughout the ordered phase. Thus this topological order survives even in the 
presence of disordered couplings. 



3 Parafermions in the chiral Z n chains 

In this section I introduce Z n -invariant spin chains and show how to write their Hamiltonian in terms 
of parafermions, Z n generalizations of fermions. 



3.1 The chiral clock/Potts model 

The most famous generalization of the two-dimensional classical Ising model is the Q-state Potts 
model. The two-state Ising spin is replaced with a "spin" with Q states, and the Hamiltonian is 
invariant under the permutation group Sq, generalizing the Ising Z2 symmetry. With this symmetry, 
the only possible nearest-neighbor interaction depends on whether the two spins are the same or 
different. While this model has the virtue of simplicity, a major disadvantage is that for Q > 4, the 
phase transition between ordered and disordered phases is not continuous |llj . 

It is thus interesting to instead consider a more general model, the "clock" model. Here each "spin" 
takes on n values, as with the Potts model, but the interactions are required only to be invariant 
under Z n symmetry. It is convenient to label the values of the spin Sj at site j by l,u),u) 2 . . . 0J n ~ , 
where u) = e 2m / n . The most general Z n -invariant coupling between two spins is then 

n-i 

-J^a m {8* j8k ) m . (12) 

m=l 

with J > by convention. To make the energy real, the coefficients must obey a* m = a n - m . The 
Z n symmetry comes from sending Sj — > cosj for all j. The S n symmetric Potts case corresponds to 
having all a m the same. 
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A traditional ferromagnetic interaction corresponds to real a m and Ja m > 0. However, allowing 
a m to be complex results in interesting behavior not possible in the Ising case. To understand this, 
consider the three-state system n = 3, and let a\ = a\ = e 1 ^ . There are three possible values of s*jSk, 
given by l,oo,co 2 = uJ. When = so that the a are real, the three values of s*SkC 1 ^ are illustrated 
in the left part of fig. [TJ For a given pair of spins, the energy is proportional to minus the real part 



O 

O 



o 



Figure 1: Illustration of the three possible values of s*s/%e 1 ^ in the Z3 model for = and 7^ 0. 

of each value, so the rightmost point in this figure Sk = Sj minimizes the energy. This indeed is a 
ferromagnetic interaction when = 0. Now take > and small, as illustrated in the right part of 
fig. [TJ The interaction still favors alignment, since the rightmost point still has the largest real part. 
However, a non-zero means that the interaction no longer is invariant if s& and Sj are interchanged, 
e.g. the spins (sj, Sk) = (1, do not give the same energy as (u), 1) do. Indeed, in the figure the two 
points on the left do not have the same real part. This means that for 7^ tt/3 times an integer, 
spatial parity symmetry in any direction is broken. For this reason, a model allowing complex a is 
referred to ELS £1 chiral clock model [T2j [13] . 

In the chiral clock model it therefore is possible to continuously interpolate between ferromagnet 
and antiferromagnet without changing / or J. The antiferromagnetic case corresponds to = 7r/3 
(times any non-zero integer), where there are (at least) two values of s*s^ that minimize the energy, 
as illustrated in the rightmost part of fig. [2] A particularly interesting case is = ±7r/6 (or vr/2), 




Figure 2: Illustration of the three possible values of s^Sfce 4 ^ for the Z3 ferromagnetic case = 0, the 
"symmetric" case = —tt/6 and the antiferromagnetic case = ±7r/3. 

"halfway" between ferromagnetic and antiferromagnetic. This is illustrated in the middle picture in 
fig. [2] As is apparent, this case is symmetric around the imaginary axis, so that for each state of 
energy E, there is one of energy — E. With the quantum Hamiltonian H to be discussed shortly, 
this means its spectrum is invariant under sending H — >• —H. The symmetric case of the chiral clock 
model shares this property with the Ising model, as discussed in the previous section. 
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3.2 The Hilbert space and the Hamiltonian 

Taking anisotropic couplings in the two-dimensional chiral clock model, one can define a quantum 
clock Hamltonian acting on a chain of L "spins" with n states each. The Hilbert space for these 
Z n -invariant chains is thus (C n )® L . The basic operators a and r generalize the Pauli matrices a z and 
a x to an n-dimensional space. Instead of anticommutation relations, the operators satisfy 



1 • 



u 

OT 



a n - x 

U TO 



-t 



r-n-1 



(13) 
(14) 



where oo = e 2m / n . Although an explicit representation of these operators is not necessary for the 
subsequent analysis, it can be useful to keep in mind the following one. Diagonalizing one of these 
two operators (say a) gives 
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(15) 



In this representation a measures the value of the spin at each site, while r shifts the spin. The 
Hamiltonian is defined in terms of operators <7j and tj acting non-trivially at site j of the chain. Each 
pair of operators (o~j, Tj) satisfies the algebra ( 13|14 ), while operators at different sites commute. 

The natural generalization of the Ising Hamiltonian to the three-state case for open boundary 
conditions is 

L L-l 



H 3 = -f 



+ r,-e 



J4> 



(16) 



There are thus three physically important parameters in ^3: f /J, 4> and 4>. The one-site term, with 
the / in front, generalizes the flip term, whereas the two-site term with J in front generalizes the 
nearest-neighbor interaction. When the phases (j) = (ft = 0, this is exactly the quantum chain version 
of the three-state Potts model. The point f = J,(ft = (ft = 0is self-dual, critical, and integrable [11!. 
It separates the ordered phase, where the Z3 symmetry is spontaneously broken, from a disordered 
phase. The spectrum is invariant under #3 — > —H3 in the symmetric case (ft = (ft = tt/6 (plus tt/3 
times any integer). 

For non- vanishing phases, the physics of the model is much more intricate (see e.g. [12\ [131 122j). 
Although this is quite interesting, only a few simple properties will be necessary for the subsequent 
analysis. The discussion in the previous subsection 3.1 of behavior of the classical model when 
the phase is changed applies here as well. Varying the phases (ft and (ft allows one to interpolate 
between the ferromagnet and antiferromagnet at constant J. Here, as opposed to the Ising chain, 
the antiferromagnetic and ferromagnetic cases are not equivalent, as apparent in fig. [2j there is only 
one way for nearest-neighbor spins to align ferromagnetically, but there are (at least) two for the 
antiferrogmagnet. Outside of the ferromagnetic and antiferromagnetic cases, the chiral Hamiltonian 
breaks spatial parity symmetry and time-reversal symmetry. 

For general n, the chiral clock Hamiltonian with open boundary conditions is 



L n—1 L—l Ti—l 

j = \ 771=1 j = l 771=1 



a, 



M°3 



(17) 
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With / and J real and non-negative by definition, the conditions on the couplings 

°*m = a n-m , Oi m = OL n -m (18) 

are necessary to make the Hamiltonian hermitian. The Hamiltonian is invariant under the 7L n symme- 
try of sending aj — > uaj for all j, i.e. "increasing" each spin. The corresponding symmetry generator, 
generalizing (—1)^, is 

which indeed satisfies (uj p ) n = 1. 

3.3 The chiral clock model in terms of parafermions 

Parafermions in an n-state spin chain are defined analogously to the fermions in the Ising/Majorana 
chain, by multiplying order and disorder parameters |15] . At each site there are two basic parafermions 
Xj and ipj, generalizing a,j and bj\ 

Xj = (lI T fcW ( 2 °) 

^ = - ( - 1)/2 X 3 r 3 = (jl r k ) a 3 r 3 . (21) 

Like a and t, these do not square to 1 and do not commute, but rather 

(Xj) n = fa)" = 1 , X) = (XjT- 1 , rf = (Tpj) n -\ (22) 
Xjipj = ui>jXj ■ (23) 

Because of the strings attached, operators at different points do not commute like Uj and do. 
Instead, 

XjXk = ^XkXj , ^j^k =^^k^j , Xjipk = uip k Xj for j < k - (24) 

The restriction j < k is necessary for these relations to make sense; only for Ising is to = . Thus 
while parafermions are a natural generalization of fermions, calculations involving them are much 
more intricate. 



The chiral clock Hamiltonian (17) can be written simply in terms of parafermions. By definition, 



t; 



Using these along with the parafermionic commutation relation (23) means that 

L n— 1 L—l n—l 

^ = -/E E « m ^ m(m_n)/2 (xj)"(^) m - J E E «-^ m(m_n)/2 (^r~ m (xj+i) m • (25) 

j=l m=l j=l m=l 

The extra factors of u look strange, but they ensure that H n is Hermitian; note e.g. that (xtyj)^ = 

ipjXj = ^Xjtp}- 

The three-state case is 

L L-l 

H 3 = -fJ2 (xJV»j ai57 + h.c. ) ~ JJ2 (V-Jxi+i aiW + h.c ) . (26) 
j'=i i=i 



9 



The symmetric case a± = a\ = e tlv ^ has a particularly simple form: 



L-l 



3=1 3=1 

Even though H n looks simple in terms of the parafermions, it does not allow for any easy computa- 



tion of its spectrum as with the Ising chain. The fundamental reason why is in the relations (24). One 



can of course Fourier transform the parafermions into momentum space as with the fermionic solution 



of the Ising chain. The catch here is that because of the j < k requirement in (24), commutation 
relations of parafermions at fixed momentum couple different momenta. Thus Fourier transforming 
the Hamiltonian does not reduce it into 4x4 blocks, the way it does for Ising. 

The chiral clock chain is however integrable for a two-parameter family of couplings ITU] . This 
integrable case is referred to as the "integrable chiral Potts chain" ( even though there is only Z n 
symmetry, not S n ); overviews can be found in |22| 123]. The two independent parameters cj> and <p are 
defined by 

n _ i<t>(2m-n) p- 1 i$(2m-n) 

sm(irm/n) sm[irm/n) 

while the flip and interaction coefficients are 

f cos(n(f)) = J cos(ra^) . (28) 

The "superintegrable" case = = n/(2n), where any value of f/J is allowed, is particularly 
interesting |14| I18j . In particular, it resembles the Ising/Majorana chain in that the terms in the 
Hamiltonian obey the Onsager algebra [JJ , the key technical observation that allowed Onsager to solve 
the Ising model without utilizing the Jordan- Wigner transformation later exploited by Kaufmann. 
Another interesting resemblance in that its spectrum is invariant under sending H — )■ —H . It turns 
out that the parafermions have a nice relation to integr ability; for example the relations ( 27p8 ) giving 



integrability also arise by demanding that there exist a parafermionic raising or lowering operator, or 
"shift mode" . I discuss this in a companion paper [21] . 

4 Parafermionic edge zero modes 

In this section I discuss an iterative procedure for finding the edge zero mode in these Z n -symmetric 
systems. This procedure works only when the interactions are chiral, i.e. the chain is not a simple 
ferromagnet or antiferromagnet. A striking result, therefore, is that only in when spatial-parity and 
time-reversal symmetries are broken does an exact edge mode exist by this construction. In the next 
section, I prove that for the appropriate couplings this procedure works to all orders, thus showing 
that this edge zero mode is exact. 

A parafermionic zero mode \l/ is an operator that commutes with the Hamiltonian but not with 
the Z n charge oj p : 

[H n ,t>]=0, u p ^ = uj^uj p . (29) 

To survive in the L — > oo limit, it needs to be normalizable: W^B = 1. The fact that it does 
not commute with w p means that it maps between different Z ra sectors. If such a mode exists, the 
spectrum therefore must be the same in all sectors. This therefore is a much stronger statement than 
simply whether a zero-energy eigenstate exists; it is a statement about the entire spectrum. 



To illustrate this, the energy spectrum in the Z3 symmetric case ((f) = (f) = tt/6 in (26)) with 
/ = J/2 is plotted in fig. [3] for open boundary conditions and just L = 4 sites. Quite obviously, 
the spectrum in each of the three Z3 symmetry sectors is very close, consistent with the existence 



10 



charge 
-1 



E 
— >■ 



X X XX X XX XXXXXXXXXXX XX X XX X X 



E = 



Figure 3: The energy levels in the symmetric Z3 case with f = J/2 and L = 4 and open boundary 
conditions, plotted so increasing energy is to the right. The overall energy scale is unimportant; 
the purpose of this plot is to illustrate the close correspondence between the spectra in different Z3 
sectors. 

of a zero mode. As L is increased, the match between the energies becomes exponentially closer. 
The alignment starts out perfectly at //J = where the edge zero mode is exact even at finite size, 
and remains excellent as f/J is increased until around .5. For higher values of f/J, the agreement 
noticeably worsens. It is hard to say whether a transition has occurred, or this is simply a finite-size 
effect. Going to large sizes means that there are far more levels, and so to distinguish them one must 
"zoom in" on a much small region of energies, thus making it difficult to distinguish whether these 
are becoming closer or not in the continuum limit. However, this plot strongly suggests that a zero 
mode survives in the symmetric case at least until f/J ~ 1/2, and very probably further. 



Figure 4: The energy levels in the Z3 ferromagnetic case with f = J/2 and L = 4, plotted in the 
same fashion as in figure [3j The energy levels for the antiferromagnetic case are given by sending 
E->- -E. 

One the other hand, in the ferromagnetic Z3 case a = a = plotted in figure [4] for / = J/2, the 
correspondence between levels in different sectors is much less strong. (The fact that the spectrum is 
identical for Z3 charges 1 and —1 has nothing to do with a zero mode, but instead is a consequence 
of the S3 permutation symmetry present only at the ferromagnetic and antiferromagnetic points.) 
Varying f/J shows that the close correspondence seems to break down for very small f/J. Although 
some of the levels in the figure are very close in value, others differ by more than 10% between sectors. 
This is a strong indication that there is no exact zero mode in the ferromagnetic case except perhaps 
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at very small values of f/J. This conclusion also applies to the antiferromagnetic case <j> = (ft = tt/3, 
since the spectrum is that of the ferromagnetic case with a sign flip. 

The existence of an exact edge zero mode is easy to show in the extreme case / = 0. In this limit, 
the spins are ordered (except at the antiferromagnetic points). Namely, each interaction energy 

n-1 

-J £ a m (a}a j+l r 

m=l 

commutes with those at other values of j, and so all their eigenvalues can be minimized. Unless the 
a m are chosen to give multiple minima (the antiferromagnet), fixing the spin at site 1 then fixes all 
the other spins in the ground state. There are then n exactly degenerate ordered ground states, with 
each of the possible eigenvalues 1, lo, . . . w n_1 of the Z n symmetry generator lo p . In fact, at / = 0, the 
entire spectrum is the same in any of the n sectors. The parafermions give a simple way of showing 
this. Similarly to the fermionic case, the operators xi and ipi do not appear in the Hamiltonian 
with open boundary conditions. Since all the terms in the Hamiltonian are of the form Xj - " 1 ^^ 
with 1 < j < L and 1 > k < L, the parafermionic commutation relations ( 23|24 ) ensure that each 
commutes with both xi and ipL- Thus 

[H n (f = 0), X i] = [H n (f = 0),i> L ]=0. 

Since 

P P Pi i p 

u Xi = u X1 W ) U ty L = U 1p L LO 

these operators change Z n sectors, xi and ipL are obviously localized at the edges, and each taken to 
the nth power gives 1. Thus xi and tfti are indeed exact edge zero modes. 

Since there is a gap, one expects that the zero mode will remain at least approximately a zero 
mode for small enough f/J. However, that is not a convincing demonstration of topological order. 
One would like to see that it remains a zero mode, exact in the L — > oo limit, in at least some range 
of f/J. This requires showing that is normalizable; recall that in the Ising/Majorana case described 



in section 2.2, the would-be edge zero mode is no longer normalizable for / > J. This is a signal that 
the phase transition is taking place. The phase diagram of the chiral clock model however is much 
more complicated [IH [T3J HH I2H and so it is not at all obvious such a simple picture holds. In 
general, without an explicit expression a precise range of f/J where the zero mode survives in the 
L — > oo limit. However, in this section I show that the more breaking of parity, the farther it survives. 
The symmetric case, "halfway" in between ferromagnet and antiferromagnet, is the case where the 
edge zero mode seems most robust. 

The major complication in demonstrating this explicitly is the fact that commuting H n with a 
parafermion does not automatically preserve "parafermion number" like it does for fermions. Namely, 



note from ( 8|9 ) that commuting the Ising Hamltonian with something linear in the dj and bj gives 
something linear in the a,- and bj. However, in the more general parafermion case, one has e.g. 

Put more prosaically, this means that there is no such thing as a free parafermion. 

The obvious way to try to construct the edge zero mode is iteratively, as with the Ising/Majorana 



chain. Namely, break the Hamiltonian with open boundary conditions (|17j) or (25) into two pieces 
H n = F + V, where F contains all the flip terms (those with an / in front), and V all the interaction 
terms (those with a J). The edge parafermion xi commutes with V, which is why it is a zero mode 
in the / = case. It does not commute with F, but instead gives 

n-1 

[H n , Xl] = [F, Xl] = fJ2 «^ m(m " n)/2 (l - u- m ){xi) n ~ m+l {ipi) m ■ (30) 
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The question is then: can this be written as the commutator of V with something else? If the answer 
is yes, then this something else can be subtracted from so that the commutator [H n , ^] yields only 
terms of order f 2 . One can then commute with F again, and attempt again to write the result as 
a commutator of V. If this iteration can be successfully repeated to all orders, then \I> is indeed an 
exact edge zero mode. 

For simplicity consider first the Z3 case, where 



[F, X i]=2if(e~ i Ui-^x\4) , 






( -e$ 






2iJ 


e-** 










J 





with a\ = a*2 = e % ^ j sin(7r/3). Since V does not involve xi> for the iteration procedure to work each 
of these two terms individually must result from a commutator with V: 

[V,X] = fa , [V,Y]= X \4 ■ 

The only terms in V possibly giving ipi on the right-hand side are those involving 4x2 and ^>ix\- 
Thus X must be some linear combination 

X = A^ 1 + Bx2 + C4x\ , 

because commuting with V maps such a linear combination onto another such linear combination: 

[V, X] = A'tpi + B' X 2 + CVjxS- This ma P from (A B ' c ) to ( A 'i B 'i c ') is linear > and so il; can be 
written in matrix form as 



(31) 



using the parametrization Si = a\ = e~ l< ^/ sin(7r/3). Answering the question is now simple: there is 
an X satisfying = ipi if this matrix can be inverted. The determinant is — 16 J 3 sin(3</>), which 

is non-zero unless / 7r/3 times an integer. Thus it is invertible and X can be found for all (f> except 
for the ferromagnetic and antiferromagnetic cases! 

The analogus matrix giving Y is minus the Hermitian conjugate, and inverting them gives 

X = 1 tyi + e 2 ^ X2 + e- 2i *u4xl) , Y = - 1 (4 + e~ 2 ^ X \ + e 2 *Wlx^) • 
4Jsm(30) " 4Jsm(30) 

Therefore the edge zero mode is to order f/J: 

*Mt = Xi ~ 2if e-*X + 2if e^ X \y + ■■■ ■ (32) 

The dimensionless expansion parameter here seems not to simply be f/J as in the fermion case, but 
rather 

/ 

2Jsin(30) ' 

The radius of convergence of this expansion is therefore (^-dependent, and its maximum value is at 
the symmetric points <j) = ir/6 (plus 7r/3 times any integer). Thus the symmetric coupling seems to 
be the most robust for the presence of topological order. 

The fact that the iterative procedure does not work at the parity-invariant ferromagnetic and 
antiferromagnetic points are approached is not a fluke of the n = 3 case; it is true for all n. The most 
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convenient method of demonstrating this is to show that the matrices generalizing that in (31) have 
zero eigenvalues. Since commuting with H n preserves 7L n charge, there are n such matrices, i.e. 

n 

[v, ^ 3 x%?] = j £ 3$ ( 33 ) 
p'=i 



for Q = 0, 1, . . . n— 1 . Each of these matrices is n x n. The commutators necessary follow from 



m- m xT + i . # xJ+J = ("- mp - u m «)tf- m+p X ]fi q ■ (34) 



The definition of the Hamiltonian in terms of parafermions (17) and the definition of (33) 
therefore gives the entries of these matrices to be 

jffi = oj m{m - n)/2 a m co- mp (l - u mQ ) , (35) 

where m = (p — p') modn. Note that the entire matrix = 0, because (tpjXj+i) an d all its powers 
commute with V . 

Finding the eigenvalues of the is straightforward because their eigenvectors are independent 
of the a m (and are the same for all Q). Labeling these eigenvectors by with r = 0, 1 . . . n — 1, 
they have elements 

=Lo qr io q{n - q)/2 , (36) 



so that the eigenvalue equation J^)v^) = \(Q' r )vp yields 

n-l 

A (0,r) = y a m cj" mr (l - uj mQ ) . (37) 



m=l 



For the three-state case, this gives for a\ = a\ = e % ^ j sin(7r/3), 

A (1 ' r) =4sin^+y(r + l)^ for n = 3, 



as can be checked directly using ((31 



Terms involving ipf (with no \2) appear for all Q ^ on the right-hand-side of (30). Thus if any 
one of the A^' r ) = for Q / 0, then the iterative procedure will not work. For the three-state case, 
this is zero only in the ferromagnetic and antiferromagnetic cases <p = ir/3 times an integer, where 
parity is conserved. For higher n, one of the a m can always be chosen to make one of the eigenvalues 
vanish. For example, for n = 4, one of A^ 1 '^ = if o?2 = iRe(Si) ± Im(Si). 

In the ferromagnetic or antiferromagnetic cases, where all a m can be taken to be real, at least 
one of the eigenvalues vanishes. This is easy to see from (37). Because a m = «m = ct n -m here, the 
coefficient of Ja m for any m / n/2 is 

co- mr {l - oo mQ ) + co mr {l - oo~ mQ ) . 

Any time Q = 2r, this vanishes for all m. Moreover, the coefficient of a n /2 is proportional to 
1 — (jj n Q/ 2 = 1 — (— 1) < 3 ) so this term vanishes for all even Q as well. Thus 

A ( 2r < r ) = for a m real . 

If parity and time-reversal symmetry are conserved, the iterative procedure for finding the zero mode 
does not work! This is reminiscent of the two-dimensional non-Abelian states found in [25J, where it 
is necessary to violate time-reversal symmetry to have a Z n topological phase. 
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5 The proof that the edge zero modes are exact 



In the preceding section I discussed edge zero modes, starting an iterative procedure for computing 
them explicitly. In this section, I show that this procedure works to all orders, provided that the 
matrices J7~W) introduced in the previous section are invertible for Q / 0. This implies that there 
exists an edge zero mode for sufficiently small f/J. I obtain a closed- form expression for only the 
"leading" piece of the zero mode, the terms that at order (/ / J) 1 ^ 1 involve xi- Nevertheless, this 
is enough information to prove that an edge zero mode exists. An added feature is that this proof 
applies when the coupling vary in space, i.e. the Hamiltonian is 

L n—l L—l n—1 

Hn = -J2T,h a ^) m -EE J ^ ^m(cr}a j+1 r ■ (38) 

j=l m=l j=l m=l 

This indicates that the edge zero mode, and presumably the topological order, is robust in the presence 
of disordered couplings. 

It is convenient to convert the problem of finding a zero mode of H n to finding a zero eigenvalue 
of an associated "Hamiltonian" %. Because the zeroth order term in Vl/ is a parafermion and the 
Hamiltonian can be expressed in terms of parafermions, the entire series can be expressed in terms of 
parafermions. Thus the edge zero mode can be written in terms of the parafermion operators, i.e. is a 
vector in the vector space V spanned by all the parafermions and their products. Since Xj = ipf = 1 
for all sites j, V is of dimension n 2L . Ordering the parafermions as Xi> V'lj Xii ^ii ■ ■ ■ j each basis 
element can be labeled as 

\P1P2 P3 ...P2L)) = X^^xf^-.-^T ■ 

where \i = . . . n 2L -l with all labels pi are interpreted mod n. The commutator of the Hamiltonian 
with any combination of parafermions can be converted into a linear operation on V, as was done in 



the last section with (31) and (33). Namely, any commutator of H n with any element of V (i.e. 



any linear combination of any products of parafermions) with the Hamiltonian can be written as 

[HnjVft] = HfivVv (39) 

for a n 2L x n 2L matrix H. An operator comprised of parafermions commuting with H (i.e. a zero 
mode) therefore corresponds to an eigenvector of T~L with zero eigenvalue. 

The explicit form of the matrix % can be found by using ( 8|9 ) in the fermionic case, while the 



general case requires (34) and 

[xT m ^T^1\ = (u- mp - ^ m )xT m+p ^T +q ■ 

The matrix % then can be written as 

% = -T - J , (40) 

L L-1 

3=1 i=i 
where %k acts non-trivially on the "state" 

v fJr =\...p,q...)) 

with p and q in the fcth and k + 1st places respectively, taking it to the states 

v v = I . . .p — m, q + m . . . )) 
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with m = 1, . . .n — 1. The corresponding non-zero matrix elements for k even are those from (35): 



(H 2j )^ = uj m{m - n)/2 a m (u~ m P - u m «) , (41) 

while those for k odd are given by replacing a m with a m : 

(«2i-i)^ = a; m ( m ~")/ 2 a m (aT"* - oj m i) . (42) 

Because the m = term does not appear above, there are no diagonal elements in %; = 0, as is 
easily checked in ( 41~|42~ ). 



The matrix Hn U is hermitian: using (42) with fx and v reversed gives 



because a n - m = a* m and oj ±mn l 2 = (— l) m . Moreover, V can be made into a Hilbert space by defining 
the inner product 

((piP2 • ■■P2L\qiq2 ■ ■ ■ q2L}) = <W<W . . . 5 P2Lq2L . 

Note that the adjoint of in this new inner product is not the conjugate of products of parafermions 
that was used to define this basis; it is just a formal definition that makes % self-adjoint. 

Thus T~L can be thought of as a "Hamiltonian" of a n 2i -dimensional quantum-mechanical system 
(recall the original system has dimension n L ). Tools such as perturbation theory familiar from 
quantum mechanics can be used to study the zero- "energy" states of the new system, i.e. the zero 
modes of the original system. In the appendix, % is written out in more standard form. This new 
"Hamiltonian" is quite reminiscent of the original Hamiltonian, seeing as it acts on a Hilbert space 
with n states per site, has a Z n symmetry, and if not for the term in parentheses, the coefficients 
would be the same. The appendix details how this "Hamiltonian" can be split into two commuting 
copies of the original (with opposite signs in front). 

I showed in the previous section that the iteration procedure works to at least to order f /J if the 



a m are such that for Q ^ is invertible, i.e. none of the eigenvalues A^' r - ) in (37) are zero. For 
the remainder of this section, assume that this condition holds. This means that %2k can be inverted 
as long as p, q, the 2k, 2k + 1 entries in the vector on which it is acting, obey (p + q) modn ^ 0. Thus 
let Q2k be this inverse, defined in general so that 

G2kH 2 k\ ...pq...)) = / H 2 kQ2k\ ■ ■ - Pq ■ ■ ■)) = <5o,(p+g)modJ ■ ■ -Pq- ■ ■)) ■ ( 43 ) 

Constructing Q 2 k explicitly is straightforward, using the eigenvectors and eigenvalues of found 
( p6p7 ). 



m 



Rephrasing the iteration procedure in terms of the action of the "Hamiltonian" T~L, the first term 
Xi in the expansion is in the new notation vq = 1 100 •••))• The only term acting on this giving a 
non- vanishing result is Hi: 

n-l 

Hv = fiHivo = fiJ2 « m w m(m - n)/2 (u,- m -l)|l-mm000...». 

m=l 

As in the last section, the question is then if this can be written as TI2 acting on some state. H2 acts 
non-trivially on the second and third slots, so its invertibility depends on the value of P2 + P3 = Q, 
which is equal to m in this state. Since m 7^ here, and by assumption H2 in the other sectors is 



invertible, a Q2 can be constructed obeying (43). Therefore to order fi/Ji, the left zero mode ^ is 



^ = ^o - ^G 2 Hiv + ... ; (44) 
Ji 
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acting with f\H\ + J\H-2 on this first-order \E' gives terms of order (/1) . The first-order expression 



for ^ is given explicitly in terms of the parafermions for n = 3 in ( 32 ) . 

One can obviously continue in this fashion. However, there is a major potential complication. In 
the first correction to ^f, there are terms involving X2- Thus H3 does not annihilate it, and so the 
next iteration 

(fin 1 + f2n 3 )n 3 g 2 'H 1 v (45) 

involves parafermions in the first four slots, i.e. those up to "02- This means that the combination 
J\H.2 + ^2^4 must now be invertible when acting on it. This not at all obviously true. In this 
particular case, by brute force one can check that it is invertible as long as H2 and individually 
are (for Q ^ 0). However, not only does brute force quickly become unwieldy, but after a few orders, 
it doesn't even work! 

The problem is familiar from degenerate perturbation theory in quantum mechanics. Let H = 
_|_ with the state \vq) be the unperturbed state, and |0 o ) any states degenerate with it 

under the unperturbed Hamiltonian H^°\ Perturbation theory starting from |?7o) works naively only 
if (O a \H^\vo) = for all a; otherwise the "energy denominator" is zero. The problem persists at 
higher orders. Denoting the state at lib. order in perturbation theory to be \vi), the naive procedure 
breaks down if the overlap (0 a \H^\vi) 7^ 0. 

In finding an zero eigenvalue state by this iteration procedure, the problem of inverting J = 
J{H.2j is essentially the same one. It has many zero eigenvalues; for example, it annihilates any 
state \p\V2Vz ■ ■ -Pl)) that has p2j + P2j+i = modn for all j = 1, ... L — 1. These zero eigenvalues 
survive even under perturbation; I show in the appendix that there remain at least n L of them. 
(Although having so many zero modes seems like good sign for having one of them be an edge zero 
mode, it does complicate the analysis!) Denote any state annihilated by J as |0 a )), and let \vi)) be 
the order (//«/) contribution to \&. Then if 

«o a |j-h»/o 

then J cannot be inverted on the appropriate subspace, i.e. there exists no v\+\ obeying 

J\v l+1 ))=F\vi)) . (46) 
This problem can be fixed in a similar fashion as in degenerate perturbation theory, by "correcting" 



the states order by order. The key point is that \vi + i)) is not uniquely fixed by the requirement (46) 



but one can add any state |0 a )) to it while still satisfying (46). Thus one can always add some linear 
combination of them to \v{)) make 

((0 6 |^(h))-^C a |0 a ))) =0 (47) 

a 

for all b. The proof of this claim is simple. Since J is Hermitian, it has a complete set of orthonormal 
eigenstates spanning the vector space V, including those with zero eigenvalue |0 a )). Since T is 
Hermitian as well, the action of J- within the set of states |0 o )) can be diagonalized, i.e. the states 
|0)) a defined so that ((0{,|.F|0 a )) = A a 5 a b. Moreover, J 7 ]^)) = ^2 a D a \0 a )) + . . . , where the states 
in the . . . have a non-zero eigenvalue of J and so are orthogonal to all the |0 a )). Thus letting 



C a = D a /A a for the states with A a 7^ 0, and C a = for the others, gives (47). Order by order, ^ 
always can be corrected so that it is an exact zero mode. 

There is a catch, however: there is no automatic guarantee that the "corrected" state still will 
be an edge state. (If there were a guarantee, then the zero mode in the ferromagnetic case could be 
"corrected" into an edge one.) In an edge state, the parafermion \i must appear in terms at least 
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order (f/J) in the expansion of ^. Nonetheless it is not difficult to show that the extra terms as 



in (47) do indeed keep an edge mode, as long as J^> is invertible for Q / 0. 

The key observation is that the iteration can be employed without any corrections to get the 
leading terms at each order. The leading terms in an edge mode at order (/ / J) 1 ^ 1 are those involving 
the parafermion xu i- e - contain only states with P21-1 7^ and pk 7^ for k > 21 — 1. I refer to such 
states as being of "length" 21 — 1. At zeroth order, the leading term is of course simply xii °f length 
1. The first-order term v\ oc Q2H1V0 contains terms of length 3, (i.e. including X2, i-e. Pz 7^ 0), as 



can be seen explicitly in the n = 3 case in (32). The second-order term V2 is found by inverting J 



acting on (45). The important thing to note is that the leading terms in (45), those of length 4, are 



not contained within Ji\V\, but rather are contained entirely within 

The reason is that all pieces of v\ have at most length 3, while %\ acts non-trivially only on the 
indices pi,P2- Thus all pieces of H1V1 still have at most length 3; only H3V1 can have pieces with 
length 4. In fact, since p^ = in all pieces of v\, all pieces of H3V1 must have P4 7^ 0. This means 
that %4 can be inverted when acting on JL3V1, because P4 7^ while = 0. Although the resulting 
state 

is only a part of V2, all the leading terms in V2 are contained within it, because the only way to get a 
state of length 5 in V2 at this order is for Ga to act on a state with length 4. Repeating this argument 
at each order means that the expression 

f\ «/l fi f\ fi fs 

^leading = VQ - — <d2H.\V§ + -— G4H3G2HW0 ~ GqHsG ' AHzG'lH.Wv + ... (48) 
J\ J\J2 J1J2J3 

contains all the leading terms in ^. 

This shows that whenever acting with %2i+i on the leading terms in vi, the resulting state has 
length 21 + 2 and that G21+2 inverts H21+2 properly; no corrections here are necessary. Corrections 
can arise only from subleading terms, i.e. those coming from acting on vi with H2j— l where j < I. 
Since terms in vi obtained from ^leading are a t most of length 21 + 1, the only states |0 a )) that can 
have 

((o a |H 2 j-ih» 7^0 , j<i, 

necessarily have length at most 21 + 1. Thus any resulting corrections necessarily have length at most 
21 + 1, and so the corrected v\ does not increase in length. The iteration procedure works! 

The only remaining question is the radius of convergence of the expansion, i.e. the largest value 
of j j J such that i s finite as L — > 00. As shown in the previous section, for n = 3, the expansion 
parameter seems to be //(2Jsin(30)). The expansion parameter for general n then seems to be 
f/(JD), where D is the smallest value of detjW = f\r=o xiQ ' r) for the 

given a m . Making / < JD 

however is not enough to demonstrate convergence, since at each order in f/J there are multiple 



terms, as is obvious from (32). The number of such terms is exponentially growing, so conceivably 
the radius of convergence is smaller; it is not possible to know exactly how much smaller without 
more explicit information about Nonetheless, for sufficiently small f/J, the series will converge. 
An interesting open question is if the convergence of this series is related to the phase diagram of the 
chain, which is much more complicated than in the Ising case |13[ 122] . 

6 Conclusion 

I have shown that there exist exact parafermionic edge zero modes in Z n -invariant spin systems for 
sufficiently small f/J, as long as the interactions are chiral. This provides a natural generalization of 
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the results for the Ising/Majorana chain, and indicate that these simple-to-define (but not to analze!) 
systems possess topological order. It would be very interesting to understand how these fit into the 
classification schemes of [91 [10] . 

A seeming major difference between the systems studied here and the Ising-Majorana chain is 
that while Id fermions exist in nature by trapping electrons in one dimension, it is hard to imagine a 
parafermion chain being formed in the same fashion. Nonetheless, Id parafermionic physics can occur 
in edges of the Read-Rezayi states for the fractional quantum Hall effect [26 i. The wave functions 
themselves are built from the correlation functions of conformal field theory [27] describing the critical 
point of the ferromagnetic chain |20j . In fact, the correlators are those of the parafermions themselves; 
in the Z2 Moore- Read case [28], these are those of Majorana fermions, given by a Pfaffian. While 
these Hall edges are purely chiral (i.e. have excitations moving in only one direction) and are gapless, 
gapped parafermionic systems resembling the ones considered here can be made by coupling two 
different Hall edges together |29[ [30] . Similar behavior also arises from coupled quantum wires j31j , 
and by including defects in multi-component Hall systems [32]. Thus it is not at all impossible for 
the topological order described here to occur in a real system. 

One particularly interesting more formal direction to explore is to understand what happens for 
closed boundary conditions. Obviously, there can be no edge modes when there are no edges, but there 
is still topological order. In the Ising/Majorana case, topological order be detected by computing the 
sign of the Pfaffian of a matrix determined by the couplings fj and Jj. This is a one-dimensional 
analog of computing the Chern number to probe topological order in two-dimensional free-fermion 
systems. No analog of this is yet known for parafermions in any dimension. Is it possible that there 
is a "Read-Rezayian" generalizing the Pfaffian here as well as in the fractional quantum Hall case? 

Another direction worthy of further exploration is to understand if similar Z n -invariant spin 
models in higher dimensions exhibit topological order. The Kitaev honeycomb model [33] in two 
dimensions arises naturally from coupling together Z2 spin chains so that the model equivalent to 
free 2d fermions in a background Z2 gauge field. Topological order results upon breaking of time- 
reversal symmetry. It is not difficult to see how to couple together Z n -invariant chains to give a model 
with a background Z n gauge field (an observation also made independently by X. Qi, and by C. Kane). 
However, since there is no such thing as a free parafermion in Id, there is not likely a notion in 2d 
either. Thus analysis of the resulting model is difficult, but perhaps the exact Id zero modes found 
here will be useful in 2d as well. At minimum, a strong hint from the results here is that it would be 
a good idea to make the interactions chiral by including these phases in the nearest-neighbor terms, 
a very natural option here. 

Finally, it would also be interesting to reverse the arrow of logic and understand what think- 
ing about topological order can say about the original spin system. A result I will discuss in a 
separate paper is that the integrable couplings ( 27|28 ) can be found simply by demanding that for 



closed boundary conditions, there exist a shift mode, generalizing the zero mode found here for open 
boundary conditions. This gives a direct and simple proof of the observations of [141 [TH] about the 
similarity in spectrum between sectors. Many other insights about the integrable spin systems using 
parafermions also appear possible. 

Acknowledgments: This work was initiated in 2010 and continued while I was visiting Microsoft 
Station Q in 2010-2011; I am grateful to Michael Freedman and the entire group for their hospi- 
tality and generosity. I am also grateful to Eduardo Fradkin, Andrei Bernevig and Taylor Hughes 
for many interesting discussions, and to Jason Alicea, David Clarke and Kirill Shtengel for discus- 
sions of their work. My research is supported by the National Science Foundation under the grant 
DMR/MPS1006549. 



19 



A Many other parafermionic zero modes 



As mentioned in section ([5]), the "Hamiltonian" 7i has at least n L zero eigenvalues. Thus the original 
Hamiltonian has at least n L zero modes comprised of parafermions, albeit with no guarantee that 
these remain normalizable as L — > oo. As also mentioned, Ti, acting on a n 2i -dimensional Hilbert 
space, is equivalent to the sum of two Hamiltonians, each equivalent (up to an overall sign) to the 
original Hamiltonian H n . The purpose of this appendix is to provide a proof of these claims. 



The matrix elements of H are given in (41) and (42). Ignoring the expression in parentheses, 
the matrix elements are exactly the same as those in H n . This suggests splitting % into two pieces 
, arising respectively from the first and second terms in the parentheses. (This is not 
the same as splitting it into flip and potential terms J- and J .) Doing this shows that not only are 
and —H^ each equivalent to the original Hamiltonian, they commute with each other. This 
allows a simple proof of the existence of exact zero modes; each corresponds to an eigenvector of 
and T-L^ with opposite eigenvalues of the two. 

Since the "Hamiltonian" T~L acts on the states \p1P2 • • •)) by simply shifting the pi, with matrix 



elements depending only on the values pi, it can be rewritten in terms of operators acting as (15). 



They thus obeying the exact same algebra ( 13|14 ) as o~j and tj do, only now acting on the new 2L 



"site" Hilbert space. To avoid confusion, I call the diagonal and shift operators s a and t a respectively. 
Using ( 41|42 ) then gives 



2L-1 n-1 



W M = - J2 «« E " m{m - n)l2 u m (t\t a+l ) m S - m , (49) 

a=l m=l 
2L-1 n-1 

u {2) = E u * E ^ m(m - n)/2 « m (tk + i) m c+i , (50) 



a=l m=l 

where U2j-i = fj and U2j = Jj and for simplicity here I set a m = a m . These two Hamiltonians 
indeed commute term by term. These new Hamiltonians are quite simple in the Ising/Majorana case. 
Since s a and t a here are the Pauli matrices a z a and a^, using a y = ia x a z gives 

2L-1 2L-1 

= " E U ^K+l > = E U a«l ■ 

a=l a=l 

The equivalence of the two Hamiltonians is obvious, coming from a canonical transformation ex- 
changing 0* and a v a at each site, sending T~L^ —T~L^. Another amusing thing to note is that if 
the canonical transformation exchanging 0* and a y a is instead done at every other sites % simply be- 
comes the famed quantum XY Hamiltonian. Splitting it into Ti^ and T~L^ then amounts in fermion 
language to splitting a free complex fermion into two Majorana theories. 

To see the equivalence between and —H^ and the original Hamiltonian, first note that all of 
the terms in T~L^ satisfy the same algebra as do the terms in the original Hamiltonian. For example, 



*a+l S a)(*a+l*a+2 s l+l) — W (*a+l*a+2 s l+l ) (4*a+l S a) ' 



just like T a {a\ L a a+ i) = u}(aia a +i)T a . The same is true for T-i^ 2 \ Moreover, the n 2i -dimensional 
Hilbert space can be split into two n^-dimensional spaces, one of which is found by acting with the 
terms in on some reference state (say 1 00 . . . ))), and the other by acting with T~L^ on the same 
reference state. Thus within each of these Hilbert spaces, n& and are equivalent to the original 
Hamiltonian H n , up to the sign in front of H^ 2 \ 

This equivalence requires HP^ 1 have the same spectrum as H n , while has the opposite spec- 
trum. Since the two can be simultaneously diagonalized, 7i = Ti^ +7i^ must have at least n L states 
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with eigenvalue zero. Thus the original problem has many zero modes comprised of parafermions! 
The harder part is now showing that one of these zero modes is a normalizable edge zero mode. This 
is the proof of section [5j 
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